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Particle  Dispersion  within  Zonal  Jets 
in  Two-Layer  Beta-Plane  Turbulence.  (August  1996) 

Jennifer  Claire  Roman,  B.  S.,  University  of  Arizona 
Chair  of  Advisory  Committee;  Dr.  R.  Lee  Panetta 

Passive  tracer  dispersion  is  studied  in  a  two-la5^er,  quasigeostrophic,  beta- 
plane  model  in  which  persistent,  steady,  zonal  jets  are  observed.  Particle  trajectories 
and  statistics  are  examined  for  barriers  to  or  mechanisms  for  mixing.  Simulations 
are  performed  for  two  different  values  of  the  planetary  vorticity  gradient,  and  an 
Eulerian  energy  timescale  and  a  lengthscale  dependent  on  the  number  of  jets  in  the 
flow  are  developed  to  enable  comparison  between  the  two  systems.  Using  this  new 
timescale  and  lengthscale,  it  is  shown  that  zonal  particle  root-mean-square  displace¬ 
ments  are  enhanced  as  j3  increases,  and  meridional  rms  displacements  in  regions  away 
from  wave-breaking  activity  are  unaffected  by  an  increase  in  the  planetary  vorticity 
gradient.  Two-particle  spreading  also  increases  with  0.  Larger  values  of  beta  result 
in  diminished  Rossby  wave-breaking  on  the  maximum  gradient  of  potential  vorticity 
associated  with  jet  cores,  so  that  tracer  dispersion  across  westerh^  jets  is  decreased. 

The  dependence  of  tracer  statistics  on  the  scales  of  waves  included  in  the 
adA^ecting  flow  is  also  examined.  Hyperviscosity,  rather  than  laminar  or  biharmonic 


diffusion,  is  used  to  parameterize  dissipation  in  the  model,  so  that  enstrophy  is  allowed 


IV 


to  cascade  to  small  scales.  As  a  result,  single-  and  two-particle  statistics  do  show  a 
dependence  on  the  number  of  small  scale  waves  included,  and  in  particular,  it  is  shown 
that  removal  of  small  scales  increases  zonal  displacements  and  inhibits  meridional 
displacements.  This  is  thought  to  be  a  result  of  the  decreased  Rossb}'^  wave-breaking 
that  accompanies  the  removal  of  small  scales  from  the  advecting  flow. 


Vll 


TABLE  OF  CONTENTS 

Page 

ABSTRACT  .  iii 

DEDICATION  .  v 

ACKNOWLEDGEMENTS  .  vi 

TABLE  OF  CONTENTS  .  vii 

LIST  OF  TABLES  .  ix 

LIST  OF  FIGURES  .  x 

CHAPTER 

I  INTRODUCTION .  1 

II  BACKGROUND  .  5 

Turbulent  Diffusion  .  5 

Potential  VUrticitj^  and  Passive  Tracers  .  9 

Rossbj'^  Wave-breaking  .  11 

III  MODEL  DESCRIPTION  . 13 

Eulerian  Flow .  13 

Lagrangian  Fields  .  16 

Statistical  Analysis .  19 

IV  SIMULATIONS  AND  DISCUSSION  .  22 

Effect  of  the  Planetar}^  Vorticity  Gradient  .  22 

Scale  Sensitivity  .  41 

V  SUMMARY  AND  CONCLUSIONS  .  69 


REFERENCES 


71 


Vlll 


APPENDIX 


Page 

73 


VITA 


77 


IX 


LIST  OF  TABLES 

TABLE  Page 

1  Parameters  for  the  two  beta  cases  .  22 

2  Parameters  used  in  scale  sensitivity  simulations .  41 

3  Lagrangian  integral  timescales  for  the  scale  sensitivity 

simulations  .  62 


X 


LIST  OF  FIGURES 

FIGURE  Page 

1  Spin- up  and  equilibration  energy  levels  .  17 

2  Zonalh?^  averaged  lower-layer  transient  eddj'^ 

zonal  flow  .  17 

3  Single-particle  displacements  for  ,5=. 25  .  24 

4  As  in  figure  3  but  for  P—A .  25 

5  Two-particle  displacements  for  ,d=.25  .  27 

6  As  in  figure  5  but  for  ^=.4  .  28 

7  Tracer  trajectories  for  /?=.25  at  lOTg  .  30 

8  Tracer  trajectories  for  /3=A  at  lOT^  .  30 

9  As  in  figure  7  but  for  20Te  .  31 

10  As  in  figure  8  but  for  20Te  .  31 

11  As  in  figure  7  but  for  30Te  .  32 

12  As  in  figure  8  but  for  SOTe  .  32 

13  Potential  vorticity  field  for  j3=.2b  at  time  step  2  .  34 

14  As  in  figure  13  but  for  time  step  3  .  35 

15  As  in  figure  13  but  for  time  step  4  .  36 

16  Potential  vorticity  field  for  (3— A  at  time  step  2  .  37 

17  As  in  figure  16  but  for  time  step  3  .  38 

18  As  in  figure  16  but  for  time  step  4  .  39 


XI 


FIGURE  Page 

19  Initial  potential  vorticity  field  for  Kco  =60  .  42 

20  As  in  figure  19  but  for  JCro=20  .  42 

21  As  in  figure  19  but  for  /Go=10  .  43 

22  As  in  figure  19  but  for  K,.„—b .  43 

23  Single-particle  displacements  for  Kco  =60  .  45 

24  As  in  figure  23  but  for  A'^o=20  .  46 

25  As  in  figure  23  but  for  Rco=10  .  47 

26  As  in  figure  23  but  for  Kco=5 .  48 

27  Two-particle  displacements  for  Kco  =60 .  50 

28  As  in  figure  27  but  for  Kco— 20  .  51 

29  As  in  figure  27  but  for  A5:o=10  .  52 

30  As  in  figure  27  but  for  /G„=5 .  53 

31  Maximum  meridional  excursions  for  Kco  =60 .  55 

32  As  in  figure  31  but  for  A'.o=20  .  56 

33  As  in  figure  31  but  for  A'c„=10  .  57 

34  As  in  figure  31  but  for  Kco=5 .  58 

35  Zonal  velocity  autocorrelations  for  the  four  cases .  60 

36  As  in  figure  35  but  for  meridional  velocities  .  61 

37  Wave-breaking  time  series  for  Kco  =60  .  63 


FIGURE 


Page 


xii 


38  As  in  figure  37  but  for  /Uco=20  .  64 

39  As  in  figure  37  but  for  Kco—IQ  .  65 

40  As  in  figure  37  but  for  Kco=^ .  66 


1 


CHAPTER  I 
INTRODUCTION 

Although  it  is  difficult  to  give  an  exact  definition  of  turbulence,  we  may 
gain  insight  into  turbulent  flow  b}'^  examining  some  of  its  characteristics;  irregular 
behavior,  intermittancy,  nonlinearity  and  strong  mixing.  Of  these,  the  most  impor¬ 
tant  feature  of  turbulence  is  its  diffusivit}'^,  or  mixing  abilit}'-  (Tennekes  and  Lumley 
1972).  Simple  dimensional  analysis  shows  that  the  rate  of  diffusion  due  to  turbulent 
motions  is  several  orders  of  magnitude  larger  than  the  molecular  diffusion  rate.  Given 
an  isotropic  (/?  =  0)  flow  with  eddy  velocity  V\  the  “eddy  diffusivity”  nj  is  defined 
by 

V'f)'  =  ktV9 

where  6  is  some  quantity  of  interest  (i.e.,  heat).  This  expression  ma}'^  be  written  as 

^T^(|U'-^1)/(|V^|) 

and  dimensional  analysis  yields 

k.t  ^  {Ue)l{eiL)  =  UL. 

The  ratio  of  the  eddy  diffusivity  to  the  molecular  diffusivity,  i/,  can  now  be  written 
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K.-rli' ^  UL/iy  =  Re 

where  Re  is  the  Reynolds  number  of  the  flow.  Thus,  the  Reynolds  number  of  a 
turbulent  flow  can  be  thought  of  as  the  ratio  of  molecular  diflTusion  to  turbulent 
diffusion,  and  in  order  for  the  Reynolds  number  to  be  large,  a  necessary  condition 
for  turbulent  flows,  u  must  be  much  less  than  kt-  In  anisotropic  turbulent  flows, 
such  as  the  ones  studied  in  this  research,  kt  is  direction-dependent,  and  its  value 
is  proportional  to  the  flux  in  a  particular  direction  divided  by  the  gradient  in  that 
direction.  Thus,  the  eddy  diffusivity  in  the  zonal  direction  could  be  very  different 
than  the  eddy  diffusivity  in  the  meridional  direction.  In  both  isotropic  and  anisotropic 
flows,  however,  turbulence  is  an  expedient  means  of  mixing  and  a  robust  mechanism 
for  transport  of  tracers  such  as  potential  vorticit}'’,  passive  contaminants  and  heat. 

The  presence  of  zonal  jets  in  a  flow  has  varying  effects  on  transport  depen¬ 
ding  upon  whether  the  tracer  is  active  or  passive  (an  active  tracer  interacts  with 
the  surrounding  flow',  wTereas  a  passive,  or  dynamical!}^  inactive,  tracer  is  merely 
carried  along  by  it) .  Eddy-forced  zonal  jets  tend  to  coincide  wdth  regions  of  increased 
meridional  transport  of  heat,  an  active  tracer.  In  the  troposphere,  the  stormtracks 
associated  wdth  jets  are  a  major  region  of  poleward  transient  eddy  heat  transport 
(Holton  1992);  the  thermal  wdnd  balance  requires  a  strong  meridional  temperature 
gradient  to  be  consistent  wdth  the  strong  vertical  shear  in  jets.  The  dispersion  of 
passive  tracers,  however,  has  been  numericall)''  (Pierrehumbert  1991a;  Bowmian  and 
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Chen  1993),  experimentally  (Sommeria  et.  al.  1989),  and  observationally  (Bowman 
1993)  shown  to  be  suppressed  in  the  meridional  direction  in  the  presence  of  zonal 
jets.  The  laminar  jets  are  sites  of  very  little  mixing  and  become  barriers  to  passive 
tracer  transport,  resulting  in  anisotropic  dispersion  with  zonal  trajectories  favored 
over  meridional  trajectories.  In  the  case  of  the  Antarctic  polar  vortex,  the  barrier  to 
mixing  plays  a  crucial  role  in  separating  the  ozone-depleted  air  over  the  pole  from  the 
ozone-rich  air  at  lower  latitudes. 

The  polar  vortex  and  other  jets  have  been  found  to  be  permeable  barriers, 
however,  due  to  waA'-e-breaking  events.  In  the  stratosphere,  Rossby  waves  break  on 
the  edge  of  the  polar  vortex,  ejecting  polar  air  into  the  midlatitudes  where  it  is  mixed 
awaj^  (McIntyre  and  Palmer  1983).  Wave-breaking  has  also  been  observed  along  jets 
in  the  upper  troposphere,  and  along  the  Gulf  Stream  (Nakamura  and  Plumb  1994). 
In  each  case,  the  transfer  of  material  across  the  jet  core  is  irreversible  (that  is  what 
defines  wave- breaking) ,  and  the  penetrability  of  the  jet  is  related  to  the  frequency  and 
type  of  wave-breaking  events. 

This  research  seeks  to  further  the  understanding  of  passiA^e  tracer  dispersion 
in  a  S3'^stem  that  combines  the  efficiency  of  transport  in  turbulent  flows,  the  general 
anisotropA^  of  transport  in  the  presence  of  jets,  and  the  episodic  nature  of  transport 
associated  with  Rossby  wave-breaking.  Efforts  are  concentrated  in  two  main  areas: 


a  description  of  zonal  and  meridional  transport,  and  the  dependence  of  transport  on 
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wave  scale. 

The  first  objective  seeks  to  qualitatively  and  quantitatively  describe  the 
dependence  of  tracer  dispersion  on  the  value  of  the  planetary  vorticity  gradient.  Due 
to  the  nondimensionalization  of  the  model,  the  model  parameter  [3  is  proportional  to 
the  inverse  of  the  vertical  shear,  i.e.,  increasing  beta  is  equivalent  to  decreasing  the 
vertical  shear,  or  decreasing  the  meridional  temperature  gradient.  In  this  research, 
two  values  of  beta  are  examined  that  equate  to  a  strong  and  moderate  horizontal 
temperature  gradient/vertical  shear.  Particle  trajectories  and  statistics  are  examined 
for  barriers  to  or  mechanisms  for  mixing,  and  efforts  focus  primarily  on  understanding 
meridional  dispersion. 

The  second  objective  describes  the  dependence  of  tracer  transport  on  the 
scales  of  waves  included  in  the  background  flow.  In  some  simulations,  only  certain 
wavenumbers  are  retained  and  experiments  are  repeated.  Removal  of  small  scale 
waves  is  equivalent  to  a  coarsening  of  grid-point  separation.  Residts  should  determine 
which  scales  dominate  particle  dispersion  for  this  S3''stem,  and  show  that  accurate 
characterization  of  tracer  transport  depends  upon  the  resolution  of  the  observational 


network. 
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CHAPTER  II 
BACKGROUND 


Turbulent  Diffusion 

Tajdor  (1921)  was  the  first  to  extract  Lagrangian  statistics  from  an  Eulerian 
turbulent  flow.  Examining  homogeneous,  isotropic  turbulence,  he  derived  a  rela¬ 
tionship  between  single-particle  displacements  and  the  time  of  particle  migration,  T, 
based  upon  velocity  autocorrelation  arguments.  TaAdor  found  that  for  short  time  in¬ 
tervals,  when  T  is  so  small  that  the  velocity  of  the  flow  changes  little  as  the  particle  is 
advected,  the  standard  deviation  of  a  tracer  from  its  initial  position  is  proportional  to 
T.  For  long  migration  times,  when  changes  in  flow  velocity  are  substantial,  tracer  dis¬ 
persion  increases  as  \/T.  Such  long-term  behavior  is  termed  “random  walk” .  Taylor’s 
(and  his  colleagues’)  work  in  the  field  of  turbulent  diffusion  was  necessarih'’  limited 
to  “  simple”  flows  since  Lagrangian  statistics  could  only  be  obtained  experimentally, 
requiring  that  measurements  be  taken  following  a  randomly  moving  fluid  particle. 
Not  until  computers  allowed  the  possibility  of  numerical^  simulated  turbulent  flows 
were  higher-order  Lagrangian  statistics  easily  obtained. 

Riley  and  Patterson  (1974)  were  the  first  to  examine  decaying,  isotropic  tur¬ 
bulence  by  numerically  integrating  the  Navier-Stokes  equation,  in  an  attempt  to  relate 
Eulerian  velocity  information  to  Lagrangian  velocity  statistics.  They  were  able  to  re¬ 


produce  Ta3dor’s  autocorrelation-based  results  for  single-particle  displacements,  but 
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more  importantly,  were  able  to  calculate  and  compare  the  velocity  autocorrelations 
for  the  tracers  and  the  flow  fleld.  Their  results  showed  that  the  Lagrangian  autocor¬ 
relation  coefficient,  R^,,  at  short  time  lags  (r  <  100 At),  is  larger  than  the  Eulerian 
autocorrelation  coefficient,  R^,  which  means  the  tracers  are  more  persistent  in  holding 
their  initial  conditions  than  is  the  flow  field.  At  long  time  lags  (r  >  1005t),  however, 
Ri  is  less  than  R^,  implying  that  the  tracers  undergo  a  greater  “memory  loss”  during 
this  stage  of  transport  than  does  the  mean  flow.  Huang  and  Leonard  (1995)  repeated 
these  results,  and  attributed  the  long-time  behavior  to  differences  in  the  energy  levels 
of  low-frequency  eddies  present  in  the  Eulerian  and  Lagrangian  frames  of  reference. 
{Re  >  Rl  implies  more  energy  is  present  in  Eulerian  low-frequency  eddies  than  in 
Lagrangian  ones.)  Huang  and  Leonard  also  examined  other  higher-order  Lagrangian 
statistics  and  found  a  self-preserving  Lagrangian  velocity  autocorrelation  coefficient 
that  allows  prediction  of  particle  dispersion  rates  (based  upon  an  eddy  length  scale, 
turbulent  velocity,  and  dissipation  rate)  in  systems  with  a  power-law  decay  of  tur¬ 
bulent  energy.  Yeung  and  Pope  (1989)  performed  several  experiments  on  station¬ 
ary,  homogeneous,  isotropic  turbulence  and  concluded  that  Lagrangian  one-particle 
statistics  are  particularly  affected  by  the  motions  of  the  large-scale  flow,  which  may 
be  sensitive  to  model  forcing  and  boundary  conditions.  (In  the  model  used  in  this 
research,  however,  effects  of  boundary  conditions  have  been  minimized).  Yeung  and 


Pope  also  showed  that  low-order  Lagrangian  statistics,  such  as  velocit}'^  autocorrela- 
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tion  coefficients  and  rms  displacements,  are  not  particularly  sensitive  to  the  Reynolds 
number  of  the  flow. 

While  the  previous  studies  examined  isotropic  turbulence,  Haidvogel  (1982, 
1985)  studied  particle  dispersion  in  numerically  integrated  homogeneous,  ,/?-plane 
turbulence.  Haidvogel  determined  that  Lagrangian  statistics  are  affected  bj'^  the  value 
of  the  planetar}'^  vorticity  gradient.  As  j3  increases,  flows  transition  from  turbulent 
to  wavelike  regimes,  and  zonally-oriented  paths  develop  in  the  streamfunction  and 
vorticity  fields.  Zonal  single-particle  dispersion  is  enhanced  with  increased  beta,  and 
dispersions  in  the  meridional  direction  are  suppressed.  This  “anomalous”  diffusion, 
or  deviation  from  random  walk  behavior,  is  also  expected  in  models  with  zonal  jets, 
with  superdiffusive  behavior  zonally,  and  subdiffusive  behavior  meridionally.  In  the 
atmosphere,  the  superdiffusivities  of  tracers  in  regions  with  zonal  jets  may  account 
for  the  long-range  transport  of  crustal  aerosols  from  their  continental  source  regions, 
in  cases  when  distances  over  which  the  aerosols  are  transported  are  much  larger 
than  those  expected  from  s\aioptic  winds  and  settling  mechanisms  (Ellis  and  Merrill 
1995).  Subdiffusive  meridional  transport  of  tracers  in  regions  of  zonal  jets  is  expected 
based  on  conservation  of  potential  vorticity  arguments  (see  below).  In  a  personal 
communication,  Hua  stated  that  her  results  (in  a  system  similar  to  the  one  used  in  this 
research)  indicated  globallj^  averaged  meridional  transport  could  best  be  described 


b}^  random  walk,  i.e.,  there  were  no  anomalies  in  transport  due  to  the  presence  of 
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jets.  These  results  may  be  due  to  statistical  averaging  of  the  broad,  weak  easterly  jet 
regions  with  the  narrow,  strong  westerly  jet  regions.  This  research  will  attempt  to 
further  quantify  inhomogeneities  in  the  meridional  diffusivity. 

Haidvogel’s  statistical  results  show  clearly  the  effect  of  beta  on  tracer  trans¬ 
port.  Lagrangian  velocit}'^  autocorrelation  coefficients  become  strongly  anisotropic 
as  the  planetary  vorticity  gradient  is  increased,  and  at  verj^  high  ,d,  the  zonal  cor¬ 
relation  does  not  reach  zero  throughout  the  period  of  integration.  The  beta-effect 
is  similarly  seen  in  rms  displacement  curves.  While  short-time  behavior  of  particle 
dispersion  on  a  /3-plane  continues  to  be  proportional  to  time  T  (as  in  the  isotropic 
case),  long-time  behavior  exhibits  anisotropj'’.  Zonal  dispersion  remains  correlated  so 
that  single-particle  displacements  grow  at  a  rate  larger  than  -v/r,  but  meridional  dis¬ 
persion  grows  at  a  slower  rate.  As  /3  increases,  the  zonal  dispacement  rate  approaches 
T,  while  the  meridional  rate  becomes  ver)'’  nearly  zero. 

Finallj^,  Haidvogel  showed  that  inclusion  of  only  the  large-scale  waves  {K  > 
10)  in  the  Eulerian  velocit}'^  fields  has  no  effect  on  the  one-particle  statistics,  implying 
that  the  rms  displacement  is  controlled  by  the  larger  eddies,  as  suggested  by  Yeung 
and  Pope.  Haidvogel’s  simulations  were  performed  using  laminar  (V^  on  vorticity) 
diffusion,  however,  so  that  most  of  the  enstrophj'-  was  dissipated  before  reaching 
small  scales.  The  model  used  in  this  research  parameterizes  the  small-scale  diffusion 
using  hj'^perviscosity  (V®  on  vorticity)  in  order  to  push  the  viscosity  to  much  smaller 
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scales  (allowing  a  true  enstrophy  cascade).  Part  of  this  research  will  be  devoted  to 
understanding  which  scales  of  motion  must  be  included  in  this  model  in  order  to 
correctly  characterize  the  tracer  dispersion. 

Potential  Vorticity  and  Passive  Tracers 

Haidvogel  also  examined  the  conservation  of  potential  vorticit}'^  following 
particle  trajectories.  Potential  vorticitj'^  (P-V-)  conservation  refers  to  the  abilit}'^  to 
account  for  all  changes  in  p.v.  along  a  tracer  trajectorj'^,  and  to  associate  the  changes 
with  effects  of  dissipation  and  forcing.  Haidvogel  found  that  tracers  experienced  un¬ 
explainable  changes  in  their  p.v.  values  shortly  after  initialization,  and  vorticity  was 
in  no  way  approximately  conserved.  His  simulation  used  a  highly  accurate  interpo¬ 
lation  scheme  for  particle  advection  to  minimize  interpolation  error,  so  discrepancies 
in  p.v.  conservation  were  attributed  to  truncation  of  the  Eulerian  flow  (used  in  the 
dealiasing  algorithm).  This  truncation  (determined  by  model  resolution)  forces  non¬ 
linear  interactions  with  wavenumbers  higher  than  Kmax  to  be  ignored,  resulting  in 
removal  of  an  important  component  of  the  Lagrangian  vorticity  balance.  This  sug¬ 
gests  that  spatial  resolution  is  important  to  the  actual  paths  the  tracers  follow,  as 
shown  by  Hua  (1994),  with  higher  resolution  resulting  in  “truer”  paths.  However, 
the  truncation  has  minimal  impact  wTen  one-particle  Lagrangian  statistics  are  com¬ 


puted,  since  these  statistics  are  controlled  by  the  large-scale  energy-containing  eddies 
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[K  <C  Kmax)-  Thus,  in  any  study  of  tracer  dispersion,  the  trajectories  of  the  particles 
may  be  in  question,  but  the  single-particle  statistics  derived  from  the  tracer  positions 
and  velocities  should  not  be. 

Pierrehumbert  (1991a,  1991b)  further  qualified  the  relationship  between  po¬ 
tential  vorticity,  its  gradients,  and  tracer  transport,  by  examining  the  effects  of  large- 
scale  two-dimensional  flows  on  the  mixing  of  tracers  and  p.v.  He  found  that  the 
characteristic  mixing  zones  and  barriers  to  transport  that  evolved  in  the  tracer  fields 
were  also  evident  in  the  potential  vorticity  fields.  Although  coherent  regions  of  p.v. 
may  resist  spreading  due  to  their  inherent  roll-up  and  rotation,  as  long  as  the}^  are  not 
so  concentrated  that  the  flow  cannot  tear  them  apart,  the}^  may  be  treated  as  passive 
tracers.  Thus,  in  contrast  to  the  traditional  view  that  potential  vorticity  gradients 
are  barriers  to  mixing,  Pierrehumbert  asserts  (1991a)  that  a  strong  p.v.  gradient  is 
a  signature  of  a  mixing  barrier  in  the  flow.  We  maj^  then  presume  that  where  strong 
zonal  gradients  in  p.v.  are  present,  meridional  dispersion  may  be  suppressed. 

In  the  model  of  zonal  jets  to  be  used  in  this  research  (Panetta  1993),  potential 
vorticity  gradients  are  strongest  in  the  core  of  westerly  jets.  This  implies  that  tracers 
should  encounter  a  barrier  to  meridional  transport  in  westerlies.  As  will  be  shown, 


however,  the  barrier  is  rendered  permeable  by  wave-breaking. 
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Rossby  Wave-breaking 

The  relationship  of  Rossby  wave-breaking  to  tracer  transport  gained  impor¬ 
tance  in  light  of  its  implications  for  ozone  chemistry.  McIntyre  and  Palmer  (1983) 
showed  that  wave-breaking  along  the  stratospheric  polar  vortex  results  in  ejection  of 
material  within  the  vortex  into  the  midlatitudes.  Using  potential  vorticity  as  their 
tracer,  Juckes  and  McIntyre  (1987)  simulated  the  process  of  wave-breaking  by  nu- 
mericalR  integrating  a  high-resolution,  one-layer  hemispheric  model  representing  the 
winter  stratosphere.  Their  results  showed  that  some  of  the  thin  filaments  of  poten¬ 
tial  vorticity,  created  by  wave-breaking,  roll  up  into  vortices,  while  some  filaments 
remain  elongated  and  act  as  if  the  p.v.  were  a  passive  tracer.  In  either  case,  high  p.v. 
air  from  inside  the  polar  vortex  was  ejected  into  lower  latitudes,  and  the  vortex  was 
eroded. 

Polvani  and  Plumb  (1992)  further  classified  wave-breaking  using  the  Contour 
Dynamics  method  on  a  single-contour  vortex.  Thej^  asserted  that  there  exists  a  forcing 
threshold  above  which  breaking  results  in  “dynamicallj^  significant”  changes  to  the 
vortex.  Below  this  threshold,  however,  breaking  results  in  weak  filaments  of  material 
being  stripped  from  the  vortex.  As  the  vortex  is  not  changed  substantially,  Polvani 
and  Plumb  label  this  sub-threshold  behavior  microbreaking.  It  may  be  surmised 
then,  that  passive  tracer  transport  across  a  jet  may  be  larger  in  cases  of  Rossbj^  wave 


breaking,  rather  than  microbreaking. 
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The  models  used  by  Juckes  and  McIntyre,  and  Polvani  and  Plumb  both 
exhibited  wave-breaking  only  on  the  equatorward  side  of  the  stratospheric  vortex,  and 
the  asymmetry  was  attributed  to  a  similar  one-sidedness  in  the  strain  rate  distribution 
(largest  on  the  equatorward  side  of  the  jet  due  to  jet  curvature).  For  this  reason, 
Nakamura  and  Plumb  (1994)  examined  wave-breaking  on  a  straight  jet,  and  found 
that  the  direction  of  breaking  depends  on  the  shear  of  the  flow  on  either  side  of  the 
jet.  In  particular,  when  the  shear  field  is  s}''mmetric  about  the  jet,  wave-breaking 
occurs  on  both  sides;  otherwise,  wave-breaking  is  stronger  on  the  side  of  the  jet  with 
the  largest  shear.  The  direction  of  wave-breaking  will,  of  course,  have  an  impact  on 
the  direction  of  tracer  dispersion,  and  subsequent  reversibility  of  the  transport. 

Recenth'’,  Bowman  (1995)  proposed  a  relationship  between  dispersion  of  par¬ 
ticles  by  wave-breaking  and  molecular  diffusion  of  the  tracers.  His  work  with  an  isen- 
tropic  stratospheric  model  showed  that  when  a  region  of  flow  that  contained  tracers 
was  repeatedly  subjected  to  wave-breaking,  the  resultant  transport  resembled  random 
walk,  or  molecular  diffusion.  He  related  the  timescale  of  transport  to  the  frequency 
of  the  wave-breaking  events,  and  showed  that  meridional  dispersion  of  tracers  in 
wave-breaking  regions  was  linear  in  time. 

It  is  clear,  then,  that  any  study  of  tracer  transport  in  a  system  that  contains 
Rossby  Avave-breaking,  must  account  for  the  frequency  of  waA'-e-breaking  events.  It  is 


unclear,  hoAvever,  hoAV  to  measure  this  quantity  in  the  model  used  in  this  research. 
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CHAPTER  III 

MODEL  DESCRIPTION 

This  research  investigates  the  dispersion  of  passive  tracers  in  and  around 
zonal  jets.  These  jets  are  simulated  using  a  two-layer,  quasigeostrophic  model  run 
under  certain  parameters.  After  the  Eulerian  flow  including  the  jets  reaches  a  sta¬ 
tistically  steady  state,  tracers  are  initialized  into  the  flow,  and  are  advected  for  50 
nondimensional  “model  daj^s.”  Tracer  positions  and  velocities  are  stored  at  each  time 
step  for  use  in  statistical  calculations. 

Eulerian  Flow 

The  Eulerian  flow  is  generated  by  a  two-layer,  quasigeostrophic,  beta-plane 
model,  such  as  that  discussed  in  Panetta  (1993)  and  Pedlosky  (1987).  The  domain  is 
doubly-periodic  in  order  to  eliminate  any  boundary  eifects  on  the  flow,  and  bottom 
topography  is  not  included.  The  two  layers  have  equal  resting  depths,  77,  and  the 
lower  layer  is  slightly  more  dense  than  the  upper  layer.  The  spatial  resolution  of  the 
model  may  be  128^  or  256^  grid  points,  superposed  on  a  domain  that  is  2'kL  Rossby 
radii  on  a  side,  where  L  equals  15  or  10,  depending  on  the  simulation.  The  large 
computational  domain  is  necessary  to  resolve  the  narrow  spectral  peaks  in  the  flow 
resulting  from  the  jets  (the  spectral  resolution  for  a  27rL  domain  is  1/T.) 


The  model  is  started  with  a  small  amount  of  energy,  introduced  by  initial- 
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izing  a  band  of  linearly  unstable  waves  whose  amplitudes  and  phases  are  centered 
spectrally  around  the  wavenumber  of  maximum  instability.  The  energ)''  sourcp  for 
“eddy”  motions,  or  deviations  from  the  mean  flow,  is  the  potential  energy  reservoir 
inherent  in  the  time-invariant,  domain-averaged  vertical  shear  in  the  zonal  wind.  In 
order  to  generate  the  zonal  jets  and  their  accompanying  “storm  tracks,”  the  model  is 
run  under  a  speciflc  set  of  parameters,  including  (a)  a  forcing  mechanism  provided  by 
the  assumed  horizontally  uniform,  unstable  vertical  shear,  (b)  a  background  gradient 
of  potential  vorticity,  and  (c)  weak  damping  that  allows  nonlinear  effects  to  dominate. 
Damping  may  be  provided  through  surface  drag  or  hyperviscosit}'-.  (In  this  respect, 
the  model  differs  from  that  of  Panetta  [1993],  which  included  biharmonic  diffusion.) 

The  potential  vorticity  for  la}^er  i,  where  f  =  1  is  the  upper  layer  and  i  =  2 
is  the  lower  layer,  may  be  written  in  terms  of  the  dimensional  streamfunction  T*  as 

Q:  =  (/o  +  +  (V^^*  +  (-i)''^i'7A') 

where  T*  =  and  the  internal  Rossby  radius  is  defined  by  A  =  [g{p2  — 

Pi)i?/(2p2/o)]^'^^-  The  evolution  equations  for  potential  vorticity  in  the  two-la}^er, 
quasigeostrophic,  beta-plane  model  can  be  written 

dQHdr  +  J(^t,  Q\)  =  (la) 

dQHdr  -f-  J(T^,  Ql)  =  -KhVHl  -  z/*  (lb) 


The  horizontal  Jacobian  is  denoted  by  J(., .),  is  the  horizontal  Laplacian  operator, 
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is  mechanical  damping  (surface  drag  or  Ekman  pumping),  and  i/*  is  small-scale 
mixing. 

Because  the  transport  of  tracers  by  the  eddy  field  is  of  interest  in  this  re¬ 
search,  the  streamfunctions  for  the  two  layers  are  decomposed  into  mean  and  eddy 
flow  components 

<I>l{x\y\r)  =  ri{x\y\r)-Uoy^  (Ic) 

(Id) 

The  streamfunctions  rp*  are  deviations  from  the  time-invariant  state  with  a  horizon¬ 
tally  uniform  vertical  shear  in  the  zonal  flow,  Uq-  It  is  assumed  the  lower  layer  is  at 
rest  (or  the  reference  frame  is  travelling  at  the  speed  of  the  lower  layer) . 

Substitution  of  (Ic)  and  (Id),  followed  bj^  nondimensionalization  (see  ap¬ 
pendix),  result  in  expressions  for  the  evolution  of  the  nondimensional  eddy  potential 
vorticity  fields 

dqi/dt  +  J(V^i,gi)  =  —dqi/dx  —  {.8  +  \)dip\ldx  —  z/V^°?/’i  (2a) 

dq2ldt  +  J(^2,92)  =  -(/?+  \)dlj)2ldx  -  (2b) 

These  equations  form  the  basis  for  the  pseudo-spectral  numerical  model.  The}^  are 
integrated  using  a  leapfrog  scheme,  with  periodic  application  of  a  leapfrog-trapezoidal 
timestep  to  eliminate  the  computational  mode.  Nonlinear  velocity  products  are  eval¬ 


uated  in  phj'^sical  space  and  then  transformed  to  spectral  space  for  integrations.  A 
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circular  truncation  is  used  and  calculations  are  dealiased  using  the  method  of  Orszag 
(1971). 

During  each  simulation,  the  model  is  allowed  to  run  several  hundreds  to 
thousands  of  model  daj'^s  until  spin-up  and  equilibration  are  complete  (see  Panetta 
1993)  and  the  zonally  aA^eraged  energy  levels  show  the  flow  to  be  statistically  steadj^ 
Figures  1  and  2  show  examples  of  the  zonally  averaged  zonal  velocity  held  and  the 
associated  energy  levels,  as  the  system  passes  through  spin-up  and  equilibration, 
reaching  a  statistically  steady  state  bj'-  the  end  of  the  period.  During  spin-up  the 
energ}'’  leA'^els  grow  rapidly,  eventually  reaching  a  state  where  levels  oscillate  about 
some  mean  value.  Correspondingly,  the  zonally  aA'^eraged  flow  undergoes  a  period  of 
jet  building  and  merging,  until  finallj'',  five  steady  jets  appear. 

Lagrangian  Fields 

Once  the  background  flow  has  reached  the  stead}''  state,  point  tracers  are 
initialized  into  each  layer  of  the  flow  field.  The  particles  are  released  with  various 
initial  configurations,  depending  on  which  part  of  the  flow  is  being  examined:  i.e.,  the 
tracers  may  be  deployed  inside  the  westerly  jets,  inside  the  easterlies,  or  on  equally- 
spaced  grid  points  throughout  the  domain. 

In  each  case,  the  particle  trajectories  are  obtained  using  the  Eulerian  eddy 
velocities  at  the  positions  of  the  tracers  (i.e.,  Lagrangian  frame  of  reference  moves 


Figure  1  Spin-up  and  equilibration  energy  levels.  Energy  magnitude  on  ordinate,  time  on  abscissa. 


Figure  2  Zonally  averaged  lower-layer  transient  eddy  zonal  flow.  Meridional  distance  on  ordinate, 
time  on  abscissa. 
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with  the  mean  background  flow): 

dxpfdt  =  v[xp,t),  p=l,mtra 

where  mtra  is  the  maximum  number  of  tracers  initialized  in  the  flow.  The  Eule- 
rian  eddy  velocity  flelds  are  interpolated  from  the  grid  points  to  the  position  of  the 
tracer  by  fitting  a  streamfunction  to  the  velocity  field  over  the  neighboring  16  (4x4) 
grid  points  and  using  a  third-order  Lagrange  interpolation  to  find  the  value  of  the 
streamfunction  (velocity)  at  the  tracer  position,  as  in  Hua  (1994).  The  evolution 
of  observed  potential  vorticity  of  each  tracer  is  also  calculated  using  the  16  point 
Lagrange  interpolation. 

Once  the  tracer  velocities  are  calculated  for  a  particular  timestep,  the  particle 
positions  are  integrated  forward  in  time  using  a  fourth-order  Runge-Kutta  scheme. 
For  successive  timesteps,  n  and  n  -}- 1,  the  position  of  the  tracer  at  n  -|-  1  is  given  by 

Xn-i-l  =  Xn  -\r  ^(^1  +  2A;2  4“  2^:3  -\-  kn) 

for 

h  =  St  v{xli,tn) 
k2  -  St  v{£ri  +  +  ^St) 

ks  =  St  v{Xn  -1-  \k2,  tn  +  \^St) 
k^  =  St  u(f^  -f  k^,  t„.  -i-  St) 
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where  St  is  the  time  step  interval.  The  tracers  are  advected  for  50  model  days  for 
each  simulation,  allowing  the  eddies  to  evolve  through  several  turnover  times. 

Statistical  Analysis 

Time  series  of  tracer  positions  and  vorticities  are  stored  during  advection, 
and  this  data  is  used  to  perform  the  Lagrangian  statistical  anabasis.  Statistics  calcu¬ 
lated  the  model  include  the  velocit}^  autocorrelation  coefficients,  root-mean-square 
single-particle  displacements  and  two-particle  relative  displacements. 

As  a  particle  moves  away  from  some  initial  position,  it  will  become  decor- 
related  from  the  conditions  existing  at  the  starting  point,  or  will  “lose”  its  memory 
of  the  initial  conditions.  A  measure  of  this  memory  loss  is  the  Lagrangian  velocity 
autocorrelation  coefficent,  which  may  be  written 


^p(^)  =  T  =  t  or  t  +  T 


where  each  integral  is  evaluated  over  all  time  lags  0  <  r  <  rnaxlag.  Brackets  () 
represent  an  ensemble  average  over  all  tracers  particles,  1  <  p  <  mtra.  The  quantity 
Qp(r)  used  in  the  autocorrelation  calculation  is  the  deviation  from  the  mean,  cal¬ 
culated  by  removing  the  average  (over  all  particles)  of  cCp  at  that  time  lag.  In  this 
research,  autocorrelations  were  calculated  separately  for  a-p  =  Up  (the  zonal  velocity 
of  the  particles)  and  Qp  =  Vp  (the  meridional  velocity)  in  each  laj^er.  Thus,  a  com- 
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parison  may  be  made  between  the  along-jet  and  cross-jet  decorrelation  times.  These 
decorrelation  times,  or  Lagrangian  integral  time  scales,  are  defined  as  the  integrals 
of  the  autocorrelation  curves  over  the  total  time  of  integration.  In  this  research,  we 
expect  zonal  velocity  autocorrelation  times  to  be  longer  than  meridional  velocit}^  au¬ 
tocorrelation  times,  since,  in  the  presence  of  zonal  jets,  it  is  easier  for  a  particle  to 
retain  its  initial  east/west  motion  than  its  north/south  motion. 

Another  useful  statistic  is  the  rate  at  which  a  particle  wanders  from  its  initial 
position.  This  single-particle  spreading  rate,  or  rms  displacement,  may  be  written 

DAf)  = 

D,(t)  =  =  ([!/,(«)  -  !/,(0)]2>'/2 

where  again,  the  angle  brackets  mean  an  ensemble  average  over  all  particles.  The 
rms  displacement  at  the  Lagrangian  integral  time  scale  may  be  used  as  an  estimate 
of  the  Lagrangian  integral  length  scale,  or  the  distance  a  particle  travels  before  it 
completely  loses  its  “memory”  of  initial  conditions.  In  a  flow  dominated  by  easterly 
and  westerly  jets,  it  is  expected  that  zonal  displacements  over  a  certain  length  of  time 
should  exceed  meridional  displacements. 

A  third  statistic  of  interest  is  the  two-particle  spreading,  or  average  relative 
separation,  defined  as 
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where  the  ensemble  average  is  taken  over  all  particle  pairs  (p,  s)  with  a  certain  initial 
separation.  This  statistic  examines  only  the  magnitude  of  separation,  not  orientation, 
and  is  controlled  by  the  flow  shear  rather  than  the  flow  velocity  (Haidvogel  1985). 

While  these  and  other  statistics  were  examined,  none  satisfactorih'^  repre¬ 
sented  the  nature  of  the  meridional  dispersion  of  tracers.  This  dispersion  will  be 
shown  to  depend  upon  the  frequenc}'^  of  wave-breaking  in  the  complex  background 
flow,  a  quantity  which  we  as  yet  have  no  means  to  measure.  The  aforementioned 
statistics  do,  however,  lend  insight  into  other  aspects  of  Lagrangian  flow. 
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CHAPTER  IV 

SIMULATIONS  AND  DISCUSSION 

Effect  of  the  Planetary  Vorticity  Gradient 

The  primary  objective  of  this  research  was  to  characterize  tracer  transport 
among  zonal  jets  on  a  beta-plane.  This  was  accomplished  bj^  examining  particle 
trajectories  and  statistics  for  the  two  simulations  listed  in  Table  1.  These  parameters 
were  chosen  to  show  the  effect  of  the  planetary  vorticity  gradient  on  tracer  dispersion. 

Since  the  value  of  beta  effects  the  energ}'’  level  of  the  flow  (energy  increases 
as  beta  decreases),  a  parameter  was  needed  that  allowed  for  relative  comparisons  of 
timescales  in  the  two  simulations.  Thus,  an  Eulerian  timescale  was  developed  by 
examining  the  eddy  kinetic  energies  of  the  two  cases.  The  eddy  kinetic  energ}''  of  a 
flow  is  defined  as 


E,  =  / E{K)dR\  0<K<  Kma. 

where  K  is  the  total  wavenumber.  The  weighted  eddy  kinetic  energy  of  a  flow  is 

=  /  A  •  E{K)dK. 


Table  1  Parameters  for  the  two  beta  cases. 


Resol. 

L 

Km 

V 

Te 

l28^ 

10 

.40 

.05 

S.AIO-^ 

1.1 

256^ 

10 

.25 

.05 

1.A10-® 

.74 
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An  Eulerian  length  scale,  Ie,  may  then  be  constructed  by  taking  the  inverse  of 
{E^/Ee).  A  corresponding  Eulerian  timescale  of  the  flow,  the  large  eddy  turnover 
time,  Tg,  may  be  defined  in  terms  of  this  length  scale  by  letting 

Te  =  \1{IeYIE,  =  IeIVe:  =  EJi^.E^)  =  yfWjE^. 

It  is  this  timescale  (for  the  upper  layer)  that  will  be  used  to  compare  the  two  cases. 
The  values  of  Tg  (in  nondimensional  days)  for  the  runs  are  also  listed  in  Table  1,  and 
it  is  immediately  evident  that  the  flow  in  the  0—A  case  evolves  more  slowly  than  the 
low-beta  case  due  to  the  decrease  in  energy  of  the  flow. 

Another  difference  between  the  two  flows,  which  must  be  taken  into  account 
for  any  accurate  comparison  to  be  made,  is  the  average  distance  between  the  center 
of  the  easterlies  (where  tracers  were  initialized)  and  the  core  of  the  westerly  jets,  Lc- 
This  distance  gives  a  characteristic  size  of  eddies  imbedded  within  the  easterly  flow. 
In  the  /?=.25  case,  Lc  is  15A  (since  there  are  two  jets);  for  p=A^  it  is  7.9A  (for  four 
jets).  The  differing  average  jet  core  separations  are  used  to  relate  dispersion  lengths 
in  the  two  simulations. 

Figures  3  and  4  show  the  upper-layer,  ensemble- average  single-particle  dis¬ 
placements  as  a  function  of  nondimensional  time  (model  days),  for  the  d=.25  and 
/?=.4  cases,  respectively.  The  average  zonal  dispersion  is  indicated  by  curve  A,  and 
curve  B  shows  meridional  dispersion.  Again,  to  correctly  compare  the  two  figures. 


the  values  for  rms  displacements  at  certain  multiples  of  Tg  were  examined.  Figure 
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3  shows  that  the  zonal  rms  displacement  for  the  /?=.25  case  is  approximately  15A, 
or  one-quarter  of  the  domain,  after  ten  turnover  times.  After  twenty  T^;,  the  zonal 
rms  displacement  for  the  same  case  is  30A,  or  almost  half  a  domain.  From  figure 
4,  it  is  seen  that  the  zonal  rms  displacement  for  the  ^=.4  case  is  10. 5A  after  ten 
turnover  times,  and  18A  after  20Te.  These  values  seem  to  indicate  that  increasing  the 
planetary  vorticity  gradient  decreases  the  average  single-particle  dispersion  in  a  flow. 
However,  when  the  zonal  displacements  are  rescaled  Lc,  the  rms  displacements 
increase  slightl}''  with  increased  beta.  After  ten  Tg,  the  values  of  zonal  displacement 
equate  to  l.OLc  for  /?=.25,  and  1.3Lc  for  8=  A.  Similarly,  after  twenty  turnover  times, 
particles  have  been  displaced  zonally  2.QLc  and  2.3Lc  in  the  low-  and  high-beta  cases, 
respectivelj^ 

Examining  the  same  two  figures,  and  concentrating  on  meridional  rms  dis¬ 
placements,  yields  different  results.  In  the  /?=.25  case,  after  ten  turnover  times  the 
average  meridional  displacement  is  7.8A.  In  the  8=  A  case,  the  meridional  displace¬ 
ment  after  lOTg  is  4.1  A.  Although  it  seems  that  particles  are  displaced  farther  in  the 
low-beta  case,  the  value  of  the  displacement  in  both  beta  cases  corresponds  to  O.SLc. 
Similarly,  after  twenty  turnover  times,  tracers  in  the  ,5=.25  case  were  displaced  lOA 
meridional^,  and  tracers  in  the  8=  A  case  were  displaced  5.2A.  These  values  corre¬ 
spond  to  approximately  O.GbL^  in  each  s5-stem.  Thus,  unless  meridional  dispersions 


are  rescaled  to  account  for  the  number  of  jets  in  the  domain  and  the  relative  energies 


29 


of  the  flows,  it  seems  that  increasing  beta  acts  to  decrease  meridional  dispersions.  The 
rescaling  by  and  Lc,  howe  ver,  indicates  that  meridional  dispersion  (in  the  absence 
of  wave-breaking)  is  not  affected  by  the  value  of  the  planetary  vorticity  gradient,  and 
zonal  dispersion  slightlj'^  increases  with  an  increase  in  beta. 

Figures  5  and  6  show  the  two-particle  displacement  curves  for  the  two  beta 
cases  as  a  function  of  model  daj'^s.  Curve  A  represents  upper-layer  spreading,  and 
curve  B  represents  the  lower  layer.  After  IOTe  in  the  upper  laj'^er  of  the  /3=.25  case, 
two  particles  have  spread,  on  average,  10. 5A,  or  0.7Lc.  After  the  same  multiple  of 
turnover  times  in  the  d=.4  case,  particles  have  spread  8.4A,  or  l.lLc-  The  values 
after  twenty  turnover  times  are  26. 5A  (l.SLc)  for  the  (5=.2b  case,  and  20A  (2.5Lc) 
for  /?=.4.  The  results  indicate  that  an  increase  in  beta  results  in  increased  particle 
spreading.  Although  two-particle  statistics  show  magnitudes  of  displacement  instead 
of  orientations,  it  is  safe  to  assume  that  the  increase  in  particle  spreading  in  the  /?=.4 
case  is  due  to  the  increase  in  average  zonal  displacement  in  that  case. 

The  next  series  of  figures  (figures  7-12)  show  actual  particle  trajectories 
superposed  on  the  zonally-averaged  position  of  the  westerty  jets.  Only  a  portion  of  the 
tracers  initialized  in  each  case  is  shown  for  clarity  of  behavior,  and  approximately  100 
timesteps  are  plotted  for  each  tracer  to  show  the  relative  movement  of  the  particles. 
Figures  7  and  8  show  the  positions  of  tracers  after  ten  turnover  times.  In  both  beta 


cases,  tracers  have  dispersed  into  the  westerlies,  but  in  the  /?=.25  case,  at  least  one 
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tracer  has  passed  through  a  jet  (lower  left  of  figure  7).  After  20Tb  (figures  9  and  10), 
more  tracers  have  crossed  the  jets  in  the  low-beta  case,  while  no  tracers  have  done  so 
in  the  (5= A  case.  Tracers  in  the  high-beta  case  have  thus  far  been  unable  to  traverse 
the  jet  core.  Finalh^,  figures  11  and  12  show  trajectories  after  thirty  turnover  times. 
Again,  in  the  low-beta  case,  some  tracers  have  been  able  to  pass  through  the  jets, 
whereas  cross-jet  dispersion  is  still  inhibited  in  the  high-beta  case.  This  sequence 
of  figures  seems  to  indicate  that,  although  meridional  dispersion  is  not  affected  by 
the  value  of  /?,  there  is  clearly  something  associated  with  the  westerly  jets  that  is  d 
dependent  and  affects  northerly/southerly  transport. 

Although  not  reproducible  here,  animations  of  the  upper-layer  potcuitial  vor- 
ticit}^  fields  clearly  show  that  Rossby  wave-breaking  is  enhanced  in  the  low-beta  case. 
Three  instantaneous  fields  from  these  animations  are  shown  in  figures  13-15  for  ,5=. 25, 
and  figures  16-18  for  0=A.  (Note;  The  low-beta  case  was  done  for  128^  resolution 
due  to  time  constraints,  so  the  quality  of  the  images  in  figures  16-18  is  degraded 
compared  to  the  256^  images.)  The  average  positions  of  the  maximum  gradient  in 
potential  vorticity  associated  with  the  jet  cores  of  the  d=.25  case  are  secur  at  lati¬ 
tudes  20  and  150  (labelled  in  grid  points).  The  Rossby  waves  in  the  low-beta  case  are 
extremely  distorted  due  to  wave-breaking,  and  the  onh''  part  of  the  maximum  p.v. 
gradient  that  appears  to  be  stable  during  this  period  of  the  integration  is  at  the  right 


side  of  the  northern  jet.  The  figures  for  the  high-beta  case,  however,  show  a  reduction 
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Figure  14  As  in  figure  13  but  for  time  step  3. 


Figure  15  As  in  figure  13  but  for  time  step  4. 
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Figure  17  As  in  figure  16  but  for  time  step  3. 
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in  wave-breaking  activity.  (The  /3=.4  case  has  four  jets,  located  at  grid-point  latitudes 
20,  55,  88,  and  116.)  Wave  amplitudes  are  much  smaller  in  this  case,  and  there  are 
several  sections  along  the  maximum  potential  vorticity  gradient  where  waves  are 
extremely  weak.  While  wave-breaking  is  still  present  when  /5=.4,  the  frequency  of 
the  breaking  and  the  scale  of  the  breaking  waves  are  significantl}'^  reduced.  Differences 
in  wave-breaking  activity  between  the  two  beta  cases  are  explainable  by  the  differing 
energies  of  the  flows.  The  total  energ}'^  in  the  /1=.25  case  is  8.6  (nondimensional 
units  of  Uo^),  and  the  0=A  case  has  a  total  energj^  of  1.5  units.  The  fact  that  the 
rescaling  of  time  in  each  beta  case,  based  upon  differences  in  eddy  energies,  resulted 
in  comparable  meridional  dispersions  indicates  that  most  of  the  eddy  energy  in  each 
system  is  contained  within  the  easterlies.  The  amount  of  total  energ}'^  in  the  flow 
does,  however,  affect  the  westerlies,  and  translates  into  increased  wave-breaking  in 
the  low-beta  case.  This  results  in  a  difference  in  tracer  transport  across  the  jets  in 
the  two  systems. 

Thus,  it  has  been  shown  that  an  increase  in  the  planetar}^  vorticity  gradient 
enhances  zonal  dispersion  and  has  no  effect  on  meridional  dispersion  in  regions  away 
from  wave-breaking,  and  increases  meridional  dispersion  near  the  westerh'-  jet  core. 
These  results  rely  upon  a  rescaling  of  the  time  and  length  scales  in  each  beta  case  to 


account  for  differences  in  energy  and  number  of  jets. 
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Scale  Sensitivity 

The  second  objective  of  this  research  was  to  examine  the  sensitivity  of  tracer 
transport  to  the  scales  of  waves  included  in  the  flow  field.  This  was  done  using 
the  system  described  in  table  2.  Contributions  from  unwanted  wavenumbers  were 
removed  from  the  background  fiow  before  advecting  the  tracers  by  multipljdng  the 
velocity  vectors  in  wavenumber  space  by  1.0  if  /C  <  Kco  and  by  0.0  if  K  >  Kco-  The 
cutoff  wavenumber,  Kco-)  was  assigned  a  different  value  for  each  of  the  four  cases  to 
be  examined:  Kco  ~  60  (i.e.,  no  waves  removed),  20,  10  and  5.  These  wavenumbers 
correspond  to  4A,  1.3A,  .6A,  and  .3A  respective^,  where  A  is  the  Rossby  radius  of 
deformation. 

Figures  19-22  show  the  effect  of  wavenumber  removal  on  the  eddy  potential 
vorticity  field  in  the  upper  la)'^er.  The  graphs  show  the  upper  laj^er  p.v.  fields  in 
which  the  tracers  were  initialized,  for  the  range  of  Kco  values.  As  the  number  of 
waves  truncated  increases  (figures  20-22),  small-scale  effects  on  the  flow  disappear, 
and  the  p.v.  fields  become  greatly  smoothed.  (Note:  for  clarity  of  the  image,  the 
contour  interval  used  in  figure  19  is  twice  that  used  in  figures  20-22.)  The  potential 
vorticity  gradient  decreases  in  the  jet  core  region,  and  the  jets  become  increasingty 


Table  2  Parameters  used  in  scale  sensitivity  simulations. 
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Figure  21  As  in  figure  19  but  for  K 
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difficult  to  pinpoint,  particularly  when  Kco  <10  (figure  22). 

In  these  scale  sensitivity  simulations,  tracers  were  initialized  in  the  easterlies 
only,  so  that  tracer  behaAuor  in  and  around  the  maximum  potential  vorticity  gradient 
associated  with  the  Avesterly  jet  core  could  be  easily  distinguished.  In  other  words,  if 
tracer  statistics  representing  meridional  dispersion  exhibit  a  decrease  of  displacement 
in  time,  we  may  conclude  that  tracers  are  experiencing  some  barrier  to  dispersion  as 
they  spread  out  from  the  easterlies. 

Figures  23-26  show  the  ensemble  average  single-particle,  or  rms,  displace¬ 
ments  as  a  function  of  time  for  the  four  Kco  values.  In  each  plot,  zonal  rms  displace¬ 
ment  is  giA^en  by  curA'-e  A,  meridional  displacement  by  curve  B.  Although  there  are 
A'-erj^  slight  differences  in  the  slopes  of  the  zonal  curA^es  in  the  four  cases,  the  east/w^est 
dispersion  in  each  maj^  be  characterized  bA'^  an  initial  period  of  super-linear  dispersion 
lasting  onty  a  feAv  timesteps,  and  a  period  of  a  few  days  in  A\ffiich  dispersion  occurs 
at  a  rate  proportional  to  T.  Tracers  spend  the  majority  of  their  time,  howeA'-er,  dis¬ 
persing  at  a  rate  slightly  less  than  linear,  though  certainly  not  approaching  the  ^/T 
behavior  expected  in  isotropic,  homogeneous  turbulence.  As  Kco  decreases,  the  long¬ 
term  behavior  of  particle  dispersion  seems  to  more  closely  resemble  linear  groAvth, 
although  the  slope  of  the  rms  displacement  curA^e  is  still  slightly  less  than  1.0.  This 
AA^ould  imph^  that  tracer  trajectories  in  the  Kco=5  case  are  more  zonally  oriented  than 


those  in  the  Kco=Q0  case.  In  fact,  after  50  nondimensional  daA'^s,  a  tracer  AA^oidd  liaA^e 
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been  displaced  an  average  of  80,  88, 100,  and  (again)  100  Rossby  radii  for  Kco—^Q-,  20, 
10  and  5,  respectively.  It  seems  that  small  scales  act  to  inhibit  zonal  single-particle 
displacements. 

When  considering  meridional  rms  displacements,  however,  the  result  is  just 
the  opposite.  Tracers  initialized  in  the  full  wavenumber  case  were  displaced  an  average 
of  20A  bj^  the  end  of  the  integration,  while  tracers  in  the  other  three  cases  were  only 
displaced  14A  in  the  same  period.  A  qualitative  comparison  of  curves  B  in  each  of 
the  cases  indicates  that  as  the  smallest-scale  waves  are  removed  from  the  flow,  the 
rate  at  which  tracers  disperse  to  the  north/south  decreases  for  the  majority  of  the 
time  of  integration.  Unlike  Haidvogel’s  SA^'stem  in  which  small  scales  had  no  effect  on 
single-particle  displacements,  it  appears  that  in  this  model,  the  small  scale  waves  act 
to  increase  meridional  dispersion,  while  hindering  zonal  displacements.  This  behavior 
suggests  a  mechanism,  which  is  enhanced  b}'’  the  presence  of  small  scales,  that  forces 
tracers  across  jets  rather  than  along  them. 

Figures  27-30  show  the  two-particle  displacements  for  the  four  cases.  Curve 
A  is  for  the  upper  la3fer  in  these  figures,  and  curve  B  is  for  the  lower  laj'^er.  The  total 
average  displacement  of  pairs  of  particles  after  the  50  daj'^s  of  integration  does  show 
a  dependence  on  the  value  of  Kco-  For  K^o  =  60,  20  ,10  and  5,  pairs  of  particles 
separate  distance  of  105,  115,  118  and  A.  This  statistic  onh^  examines  the  magnitude 


of  particle  spreading,  not  the  orientation,  so  nothing  can  be  said  about  the  direction 
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of  particle  displacements.  Since  the  particle  spreading  rate  is  controlled  the  shear 
of  the  flow  rather  than  the  flow  velocity,  this  suggests  that  inclusion  of  contributions 
to  the  shear  b}?^  small  scale  waves  are  necessary  to  accurately  portraj'^  the  nature  of 
two-particle  spreading.  Thus,  while  Haidvogel  found  that  smoothing  of  the  velocity 
fields  did  not  affect  two-particle  statistics,  results  from  this  model  indicate  that  small 
scales  do  contribute  substantiallj'^  to  the  flow  shear. 

The  decrease  in  meridional  dispersion,  seen  in  the  rms  displacements,  as  the 
cutoff  wavenumber  dropped  from  greater  than  or  equal  to  10  to  less  than  10  is  also 
reflected  in  the  maximum  north/south  excursion  experienced  by  any  particle  in  each 
simulation.  Figures  31-34  show  these  maximum  excursions  for  the  first  30  days  of  each 
case,  with  solid  lines  representing  the  greatest  northward  distance  travelled  by  any 
particle,  and  the  dashed  lines  representing  the  same  for  southward  movement.  For 
the  cases  with  Kco  >  10,  there  is  no  indication  of  any  barrier  to  meridional  movement, 
although  the  rate  of  increase  of  maximum  excursion  decreases  slightty  as  more  waves 
are  removed  from  the  flow.  The  difference  in  maximum  excursion  between  the  /<’co=60 
and  JTco=20  case  is  approximately  2A,  or  3Ax,  where  Ax  is  the  space  between  grid 
points.  Thus,  removal  of  wavenumbers  larger  than  10  has  a  slight  effect  on  meridional 
dispersion.  In  the  Kco=5  simulation,  however,  there  is  clearl}^  some  mechanism  that 
acts  to  minimize  meridional  transport,  so  that  by  the  end  of  the  integration,  tracers 
in  this  case  have  dispersed  only  half  the  distance  of  tracers  in  the  Kco=Q0  case. 
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Figure  31  Maximum  meridional  excursions  for  /1l’co=60.  Maximum  displacement  in  A  on  ordinate, 
timesteps  on  abscissa.  Solid  line  indicates  northerly  excursions,  dashed  line  indicates  southerly 
excursions. 
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Figure  32  As  in  figure  31  but  for  Kco  =  20. 
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The  importance  of  small  scale  waves  to  tracer  transport  statistics  is  also  seen 
in  the  velocity  autocorrelation  calculations  for  the  four  cases.  Figure  35  shows  the  cor¬ 
relation  of  particles’  zonal  velocities  (in  the  upper  layer)  with  their  initial  conditions, 
thru  45  days  of  integration.  As  the  previous  figures  in  this  section  have  shown,  results 
for  Kco—^  are  dramatically  different  than  for  cases  with  higher  cutoff  wavenumbers. 
Velocities  of  particles  in  the  low  cutoff  case  are  well  correlated  with  initial  conditions, 
indicating  that  the  majority  of  particles  spend  their  time  in  the  easterlies.  The  other 
three  curves  in  the  figure  show  the  zonal  velocit}'^  autocorrelations  for  the  Kco  >10 
cases.  The  curves  are  very  different,  and  in  some  instances,  out  of  phase  with  each 
other  (cutoff  wavenumbers  60  and  20  at  day  7,  for  example).  Particle  velocities  for 
these  wavenumbers  rapidly  decorrelate  from  initial  conditions,  with  the  full  simula¬ 
tion  decorrelating  the  fastest.  In  fact,  for  the  case  with  all  scales  included  (Ffco=60), 
the  average  position  of  particles  initialized  in  the  easterlies  is  often  in  the  westerlies. 
This  situation  never  occurs  with  the  or  5  simulations,  and  onty  occurs  for  the 

Kco— 20  case  after  30  days  of  integration.  This  suggest  that  the  small  smales  play  a 
key  role  in  the  mixing  of  particles. 

Figure  36  shows  the  same  statistic  for  the  meridional  velocities  of  the  par¬ 
ticles.  Again,  all  of  the  curves  are  different,  indicating  that  small  scales  play  an  im¬ 
portant  role  in  tracer  transport.  Unlike  the  zonal  autocorrelations,  the  meridional 


velocity  autocorrelation  curves  are  centered  around  zero,  and  each  zero  crossing 
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Figure  35  Zonal  velocity  autocorrelations  for  the  four  cases.  Correlation  coefficient  on  ordinate, 
nondimensional  time  on  abscissa. 
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Figure  36  As  in  figure  35  but  for  meridional  velocities. 
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reflects  a  reversal  in  the  average  direction  of  movement  of  particles  compared  to 
their  initial  conditions.  Again,  there  are  instances  in  the  meridional  flow,  when  the 
movement  of  particles  in  two  cases  is  of  opposite  sign  (cutofiF  wavenumbers  60  and  5 
at  days  7  and  12,  for  example). 

Table  3  lists  the  Lagrangian  integral  timescales  obtained  from  the  autocor¬ 
relations  coefficient  curves  for  the  45  daj'^  integration  period.  Again,  the  results  reflect 
the  dependence  of  tracer  statistics  on  cutoff  wavenumber. 

Thus,  while  Haidvogel  (1985)  claimed  that  only  wavenumbers  less  than  or 
equal  to  10  needed  to  be  included  in  the  flow  to  accurately  represent  tracer  transport 
statistics,  results  for  scale  sensitivity  in  the  model  used  in  this  research  indicate  that 
for  accuracy  in  some  statistics,  small  scale  effects  must  be  included.  A  connection  be¬ 
tween  Rossby  wave-breaking  and  particle  dispersion  was  made  in  the  previous  section 
of  these  results,  and  this  relationship  maj^  explain  the  effect  of  scale  selection  on  tracer 
transport  in  these  four  cases.  Figure  37  is  a  time  series  (at  half-daj''  increments)  of  a 
portion  of  the  upper-layer  potential  vorticity  field  for  the  full  wavenumber  simulation. 
There  are  three  examples  of  waAi^e-breaking  in  this  segment  of  the  domain.  The  large 
trough  visible  in  the  left  side  of  the  upper  panel  folds  up  on  itself  and  releases  a 

Table  3  Lagrangian  integral  timescales  for  the  scale  sensitivity  simulations. 
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filament  of  positive  potential  vorticity  (solid  lines)  to  the  southwest.  The  “M”  shaped 
feature  to  the  left  of  center  in  the  upper  panel  compresses  and  releases  a  vortex  south 
of  the  jet.  Finally,  a  wave  on  the  right  side  of  the  upper  panel  has  alreadj^  collapsed, 
and  is  also  ejecting  a  vortex  into  the  flow  south  of  the  jet.  Figure  38  is  the  same 
time  series  with  wavenumbers  larger  than  20  removed.  The  three  examples  of  wave¬ 
breaking  are  still  clearty  seen,  but  there  is  a  noticeable  “smoothing”  of  the  potential 
vorticity  contours.  (Again,  the  contour  interval  in  figure  37  is  twice  as  large  as  the 
interval  used  in  figures  38-40  for  clarity.)  Some  features  are  lost  entirely  in  the  Kco=20 
case,  such  as  the  filament  released  from  the  southern  side  of  the  forming  vortex  in  the 
middle  of  the  second  panel  (see  full  wavenumber  case).  As  the  cutoff  wavenumber 
is  further  reduced,  the  wave-breaking  is  reduced.  Figure  39  shows  the  time  series 
for  the  ATco^  10  case.  Here,  the  vortex  on  the  right  side  of  the  panels  never  forms, 
eliminating  one  method  of  transporting  tracers  southward  across  the  jet.  Similar!}', 
no  clear  filament  is  formed  from  the  breaking  wave  on  the  left  side  of  the  panels, 
further  reducing  the  means  for  meridional  dispersion.  The  middle  vortex  evident  in 
the  full  wavenumber  case  also  never  forms  in  the  Kco=^^  case,  although  the  potential 
vorticity  gradient  in  that  region  does  orient  itself  in  a  north/south  direction,  allowing 
meridional  transport  of  any  tracers  on  either  side  of  the  gradient.  Finally,  figure  40 
shows  the  time  series  for  the  Kco—^  case.  In  these  panels,  there  are  no  filaments 
or  vortices  formed,  and  wave-breaking  takes  on  the  form  of  regions  of  positive  or 
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negative  potential  vorticity  becoming  disconnected  from  each  other.  The  center  of 
the  bottom  panel  shows  a  wave  reforming,  i.e.,  two  regions  of  positive  p.v.  are  being 
reconnected.  There  is  certainly  no  “crashing”  of  waves  with  subsequent  ejection  of 
material,  however. 

It  seems,  then,  that  the  importance  of  a  particular  wavenumber  to  tracer 
transport  may  depend  upon  the  wavenumber’s  contribution  to  Rossby  wave-breaking 
along  the  maximum  potential  vorticity  gradient  at  the  westerly  jet  core.  Including 
only  the  large  waves  in  the  flow  eliminates  much  of  the  breaking  and  results  in  de¬ 
creased  meridional  dispersions.  It  is  also  not  clear  that  including  only  the  large  waves 
gives  accurate  results  for  tracer  statistics.  It  appears,  then,  that  the  model  used  in 


this  research  is  sensitive  to  the  scale  of  waves  included  in  the  flow. 
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CHAPTER  V 

SUMMARY  AND  CONCLUSIONS 

This  research  characterizes  tracer  dispersion  in  a  system  that  combines  the 
efficiencj'’  of  turbulent  mixing,  the  anisotropy  of  zonal  jets,  and  the  unpredictability 
of  Rossby  wave-breaking  events. 

Results  of  simulations  with  two  different  values  of  the  planetarj'^  vorticity 
gradient  showed  the  need  to  develop  a  new  time-  and  lengthscale  to  enable  compar¬ 
ison  of  the  disparate  regimes.  An  Eulerian  timescale  was  devised,  based  upon  the 
values  of  the  eddy  kinetic  energy  and  weighted  eddy  kinetic  energy  of  the  flow.  The 
new  lengthscale  took  into  account  the  difference  in  average  distance  from  easterly 
to  westerly  jet  cores  in  the  two  beta  cases.  Anatysis  of  tracer  dispersion  statistics 
and  trajectories  using  these  new  scales  indicates  that  an  increase  in  beta  results  in 
increased  zonal  transport,  while  meridional  transport  away  from  the  westerly  jet  core 
is  comparable  in  the  two  cases.  Increasing  p  also  results  in  a  increase  in  particle 
spreading.  Similarly,  since  the  energy  of  the  flow  decreases  as  p  is  increased,  Rossby 
wave-breaking  is  diminished  in  the  high-beta  case,  resulting  in  diminished  tracer 
transport  across  jets. 

The  scale  sensitivity  experiments  confirm  that  wave-breaking  plays  a  key 
role  in  tracer  dispersion.  As  the  number  of  waves  removed  from  the  flow  increases. 


the  maximum  potential  vorticity  gradient  associated  with  the  jet  core  is  smoothed 
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and  small-scale  features  of  the  flow,  such  as  filaments  and  vortices,  no  longer  ap¬ 
pear.  Wave-breaking  diminishes  and  the  jets  become  less  permeable.  As  a  result, 
zonal  single-particle  displacements  are  enhanced  and  meridional  rms  displacements 
decrease.  Zonal  velocities  of  particles  remain  correlated  with  initial  conditions  much 
longer  when  small  scales  are  removed,  and  on  average,  the  zonal  and  meridional 
Lagrangian  integral  time  scales  increase  as  the  cutoff  wavenumber  decreases.  The 
removal  of  small  scales  also  has  an  effect  on  the  two-particle  statistics,  and  results  in 
an  increase  in  particle  spreading  as  waves  are  truncated  from  the  flow.  Thus,  the  scale 
sensitivity  simulations  confirm  that,  for  systems  where  enstrophy  is  dissipated  by  h}^- 
perviscosity,  all  scales  of  waves  must  be  included  in  the  flow  in  order  to  accurately 
characterize  tracer  dispersion. 
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APPENDIX 


The  dimensional  evolution  equations  for  potential  vorticity  are 

+ (i) 

dQHdt*  +  J(T^,  Ql)  =  (2) 

Working  with  the  left  hand  side  of  equation  (1) 
dQHdV  +  Ql)  =  Jr(V2^r*  _  #7A2)+ 

(Mt/5x*)(5gi/5r)  -  {d^l/dy*){dQUdx^)  (la) 

=  dql/dr  +  diplldx*{0*  +  -  T7A2))- 

{drjdy-  -  Uo)i,{V^'^l  -  4>*/A^))  (lb) 

=  dq^/dr  +  0*d^pt/dx-  +  (d^;/dx*)^(v^^*,)-(d^^/dx^)(i/2A^)^(^^- 
Uoy*  -  ^2)  -  +  {dipildy*){l/2\^)-^{ipl- 

Uoy*  -  ^2)  +  -  CA)(i/2A2)^(^*  -Uoy-  0;)  (Ic) 

=  dql/dr  +  P*drjdx*  +  idijl/dx*)^{VH^  -  {dryldx-)^{riX^)  + 
{dri/dx*){Uol2X^)  -  {dxl/lldy*)i,iy^xiyi)  +  {drj dy*)£:{r  ! X^)+ 
Uoi.{^^r,)-Uoi,{rix^)  (Id) 

=  dq*Jdr  +  {drjdx*)^{V^ry  -  r/A")  -  {drjdy*)^{V^^ri- 

rix^)  +  {drjdx^p*  +  6^/2A2)  +  U^{dqlldx*)  (le) 
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So, 

dqlldt*  +  J{ru  ql)  +  {drildx*){P*  +  Uo/2\‘^)  +  Uoidq*Jdx*)  =  RHS  (If) 

Or, 

dql/dt*  +  J{r„ql)  = 

-Uo{dq*Jdx*)  -  {drxldx*){p*  +  C/o/2A2)  -  (3) 

Working  with  the  left  hand  side  of  equation  (2) 
dQ’Jdf  +  Ql)  =  -  *7^2)+ 

{d<Slldx-){dQlldy)  -  (aVJdy-){dQlldx-)  (2a) 

=  dqiidf  +  a<pi/dx-{D-  +  -  'i-VA^))- 

{di>l/dy-)^{VHl-^-l\^))  (2b) 

=  aqiidf+a-ariidx- + {drjoxQ^iv^^ii) + (dr^idx-)(\i2\^)^(r,- 

Uoy’-r,)  -  (artiay-)i:(v^n)  -  (ar2iay-)(\i2\^)^(r,- 

Uoy"  -  i>l)  (2c) 

=  dqiidf  +  0-dr2/ax-  4  (a<i,iiax-)^(v^r2)  +  (si>iiax’)^(ri\^)- 

{a^j/ai-)(t/„/2A2)  -  (ai>iiay-)^(v^r2)  -  {ar2/ay-)i:(’f-/>'‘‘)  (2d) 

=  aqi/at-  4  (ar2/ax-)^(v^r2  +  r/x^)  -  {ar2iay)-^(v‘‘r2^ 

rix^)  4  {ar2/ax-){0-  -  c/„/2A7 


(2e) 
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So, 

dqlldt*  +  Jir^.ql)  +  {dr2ldx*){p*  -  C/o/2A2)  =  RHS 
Or, 

dqlldt*  +  J{r^,ql)  = 

-{dr2/dx*){p*  -  t/o/2A2)  - 
Thus,  the  dimensional  equations  may  also  be  written  as 

dqi/dr  +  j{ruql)^ 

-Uoidqlldx*)  -  +  C/o/2A2)  - 

dqildt*  +  j{r2.qi)  = 

-{dr2/dx^){P*  -  To/2A2)  - 
Using  the  following  nondimensionalization  parameters 

Q-  =  QiUol\  'h*  =  ^.AUo, 

13*  =  /?t/o/A2,  =  kmU^IX, 


equation  (3)  may  be  written  as  (with  boxes  indicating  dimensions) 


t/o/A 

A/t/„ 

dqildt+ 

\Un  Uo  ! ^ 

A  A 

Jii’uqi)  =  - 

{dq.ldx)- 


Ujl 

A2 

(^+i) 

U„\ 

A 

{dij^ildx)- 

UoX^ 

u 

Uo>> 

Alo 

ui 

A2 


dqi  /  5t+ 


J{vuqi) 


{dqildx)~ 


{(3  +  \){d^ildx)- 


(2f) 


(4) 


(3) 

(4) 


(3b) 


(3c) 


and  equation  (4)  may  be  written  as  (with  boxes  indicating  dimensions) 


5^2/5^+ 


J{'fp2,q2)  =  - 


{p  -  l){dxlj2ldx)- 


Dividing  (3c)  and  (4c)  ^  gives  the  nondimensional  equations: 

dqi/dt  +  J{yji,qi)  -  —dqi/dx—  {(]  +  ^)dipifdx  - 

dq2l dt-\r  J {xp2, 92)  =  -(/d  +  \)dij)2ldx  - 


I 
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